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ABSTRACT 

We present new high resolution spectroscopy from which we derive abundances and 
radial velocities for stars in the field of the open cluster Tombaugh 2, which has been 
suggested to be one of a group of clusters previously identified with the Galactic Anti- 
center Stellar Structure (also known as the Monoceros stream). Using VLT/FLAMES 
with the UVES and GIRAFFE spectrographs, we find a radial velocity (RV) of 
(Vr) — 121 ± 0.4 km s~^ using eighteen Tombaugh 2 cluster stars; this is in agreement 
with previous studies, but at higher precision. We also make the first measurement of 
Tombaugh 2's velocity dispersion, which is aint = 1-8 ± 0.3 km s^^. Our abundance 
analysis of RV-selected members finds that Tombaugh 2 is more metal-rich than pre- 
vious studies have found; moreover, unlike the previous work, our larger sample also 
reveals that stars with the velocity of the cluster show a relatively large spread in 
chemical properties (e.g., A[Fe/H] > 0.2). This is the first time a possible abundance 
spread has been observed in an open cluster, though this is one of several possible ex- 
planations for our observations. While there is an apparent trend of [a/Fe] with [Fe/H], 
the distribution of abundances of these "RV cluster members" also may hint at a pos- 
sible division into two primary groups with different mean chemical characteristics — 
namely (([Fe/H]), ([Ti/Fe])) ~ (-0.06, +0.02) and (-0.28, -1-0.36). Isochrone fitting to 
the colour-magnitude distribution of apparent Tombaugh 2 members yields an age of 
2.0 Gyr, E{B-V) = 0.3, and (m-M)o = 14.5 or d = 7.9 kpc for both populations — 
parameters that are within the range of previous findings. Based on position and kine- 
matics Tombaugh 2 is a likely member of the GASS/Monoceros stream, which makes 
Tombaugh 2 the second star cluster within the originally proposed GASS/Monoceros 
family after NGC2808 to show some evidence for internal population dispersions. How- 
ever, we explore other possible explanations for the observed spread in abundances 
and two possible sub-populations, with the most likely explanation being that the 
metal-poor ([Fe/H] = —0.28), more centrally-concentrated population being the true 
Tombaugh 2 clusters stars and the metal-rich ([Fe/H] = —0.06) population being an 
overlapping, and kinematically associated, but "cold" (cry < 2 km s^^) stellar stream 
at Rgc > 15 kpc. 

Key words: Galaxy: open clusters and associations - Galaxy: fundamental param- 
eters - Galaxy: structure - Galaxy: disc - Galaxy: open clusters and associations: 
individual (Tombaugh 2) 
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1 INTRODUCTION 

Tombaugh 2 (To2; at Galactic coordinates [l,h] — 
[232.8,-6.9]°) is an old, distant open cluster in the outer 
Galactic disc. We originally targeted this cluster as part of 
an ongoing program (e.g., .Carraro et al. 2007) to expl ore in 
depth a set of star clusters proposed by iFrinchabov et al.l 
l|2004l — on the basis of their position and kinematics 
— to be possibly associated with the Monoceros stream 
(Mon; Newberg et al. 2002, Ibata et al. 2003, Yanny et al. 
2003), also known as the Galactic anticenter stellar struc- 
ture (GASS; Rocha-Pinto et al 2003, Crane et al. 2003). 
GASS/Mon was discovered as an overdensity of stars near 
the Galactic plane that seems to wrap around the outer 
parts of the Galactic disc and is frequently explained as 
tidal debris from the disruption of a dwarf galaxy on a low 
inclination orbit (e.g., Crane et al. 2003, Yanny et al. 2003, 
Penarrubia et al. 2005). To2 has also been associated with 
the reputed dwarf galaxy in Canis Major (CMa; Bellazz- 
ini et al. 2004), which is suggested to be a progenitor of 
the Mon/GASS structure (Martin et al. 2004a). However, 
Rocha-Pinto et al. (2006) argue that a more likely progen- 
itor of Mon/GASS may lie somewhere in the region of the 
sky covered by the former constellation Argo Navis, and that 
the overdensity reported near CMa is correlated to a win- 
dow in the dust extinction there; in this case, To2 would be 
farther displaced from the putative Mon/GASS progenitor. 
In any of these circumstances To2 is interesting as a poten- 
tial member of a star cluster system accreted from another 
galaxy by the Milky Way. 

Subsequently, our spectroscopic analysis of this cluster 
reveals To 2 to be interesting in its own right as the first 
known open cluster exhibiting possible evidence for multi- 
ple populations or an internal population spread, evidence 
we present and explore here No well studied open cluster 
has shown a spread in abundances of major elements (e.g., 
iRandich et alTl2006l . study of M67). The possible associa- 
tion of To2 to Mon is a characteristic shared with another 
unusual star cluster showing multiple populations, the glob- 
ular cluster NGC 2808 (Crane et al. 2003, Martin et al. 
2004a) Q But whereas a few other globular clusters with 
multiple populations are now known (see §6.2), to date no 
multi-population open clusters have been reported. Never- 
theless, having two unusual, multi-population clusters po- 
tentially associated with the GASS/Mon/CMA/Argo struc- 
tures is intriguing and may point to a mechanism for their 
creation ... if we can first establish that association and de- 
termine what GASS/Mon/CMA/Argo is. 

The nature and/or reality of the proposed Mon "tidal 
debris stream" and the CMa overdensity have been called 

expressed in this material are those of the author (s) and do not 
necessarily reflect the views of the National Science Foundation. 
X Based on observations collected at the European Southern Ob- 
servatory, Chile; Proposal 076.B-0263 

^ It should be noted that Casetti-Dinescu et al. (2007) have ex- 
cluded NGC 2808 from being a part of the Canis Major/Mon 
structures based on the model of Penarrubia et al. (2005). How- 
ever if Monoceros/GASS and Canis Major are not related or if 
NGC 2808 is like uj Cen (as suggested by Casetti-Dinescu et al. 
2007) and has its own tidal stream, the derived NGC 2808 orbit 
remains consistent with the tilted spatial configuration found in 
outer disc star clusters (Frinchaboy et al. 2004). 



into question and are currently a matter of great debate. 
Ibata et al. (2003) originally proposed that Mon could be 
related to warps in the outer disc, while Momany et al. 
(2004, 2006) have argued that much of the observed stel- 
lar overdensity associated with Mon — and particularly all 
of that associated with CMa (at / ~ 240°) — is due to 
the warping and flaring of the Galactic disc, and that no 
"extra-Galactic" component is needed to account for the 
apparent overdensities in the third quadrant. This conclu- 
sion has been contested by Martin et al. (2004b) on the basis 
of radial velocities (but cf. Momany et al. 2006), while the 
discovery of "blue plume" stars in this part of the sky has 
been used to argue further for the presence of a dwarf galaxy 
nucleus in CMa (Bellazzini et al. 2004, Martinez-Delgado et 
al. 2005, Dinescu et al. 2005, Butler et al. 2007). However, 
these young stars have also been more prosaically attributed 
to the presence of spiral arm structure (Carraro et al. 2005, 
Moitinho et al. 2006). And Grillmair (2006) has quite clearly 
identified tidal streams across the Galactic anticenter region 
with relatively low inclination (35°) and at a similar distance 
and main sequence turnoff colour as Mon, though he con- 
cludes on the basis of the best fit orbits that these streams 
have nothing to do with Mon or CMa. 

What is clear from all of the debate, and a commonly 
expressed sentiment by all sides, is that further study, partic- 
ularly detailed spectroscopic analysis, is needed before the 
complicated nature of the outer disc, its warp, fiare, spi- 
ral arm structure and possibly co-located tidal debris and 
halo substructure, can be confidently disentangled. Thus, a 
number of studies (e.g., lYong et al]|2005l : ICarnev et al.ll2005l : 
IFrinchabov et a r 2006) are being conducted to explore fur- 
ther the kinematical and chemical properties of Monoceros 
and the outer Galactic disc. Our own venture in this regard 
starts by focusing on high resolution spectroscopic investi- 
gation of old, open star clusters of the outer disc, includ- 
ing those hypothesized t o be p arts of the debated overden- 
sities. In ICarraro et all (|2007l ) we present high resolution 
spectroscopy of the first five outer disc open clusters ob- 
served in our survey (Ruprecht 4, Ruprecht 7, Berkeley 25, 
Berkeley 73 and Berkeley 75) and derive detailed abundance 
analyses and kinematics for these systems. Here, we discuss 
separately the case of the open cluster To2, which has sur- 
faced as a potentially unusual system among the clusters of 
the outer disc. In this paper we obtain new estimates of the 
radial velocities (§3) and abundances ([Fe/H] and [Ti i/Fe]; 
§4) of stars in the To2 field. We discuss the unusual abun- 
dance results for To2 in §5, and address the possible origins 
of the observed mixed chemistry, its implications for To2 
and the outer disc, and the issue of the possible connection 
of To2 to the Monoceros/GASS structure in §6. 



2 OBSERVATIONS AND DATA REDUCTION 

Our spectroscopic data of To2 come from the ESO- 
VLT (proposal 076.B-0263), a nd were collected with 
the FLAMESQVLT-HGI RAFFE j Pasquini et al.ll2002l ) and 
FLAMES@VLT-HUVES (iDekker e t al. 200(J) spectrographs 
on the nights of UT l-Dec-2005 and 5-Dec-2005. The sky 
was clear, and the typical seeing was 1.0 arcsec. FLAMES 
allows the use of the GIRAFFE spectrograph to obtain RVs 
for up to 132 stars at one time, while simultaneously using 
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Figure 1. Radial velocity distribution of all stars observed with the VLT. The cluster is clearly visible as the peak near ~ 121 km s~^. 



the UVES multi-fiber mode, which allows observing up to 
eight additional sources. For UVES, we used the 580nm set- 
up {R = 40, 000 in the 4750-6800 A range). The eight fibers 
were placed on the brighter probable members of the cluster 
(selected according to their position in the colour-magnitude 
diagram; CMD), while FLAMES-(-GIRAFFE was used to 
target 35 additional probable members, with the remaining 
FLAMES -f GIRAFFE fibers samphng the sky. In Table[T]we 
report properties of the observed stars: the first and second 
column give their ID, the third and fourth the coordinates, 
the fifth and sixth the V magnitude and V — I colour, the 
seventh, eighth, and ninth columns the heliocentric radial 
velocity with the error and the membership classification 
(see below), while the last three columns give the adopted 
atmospheric parameters used for abundance measurements 
for the member stars (see §3.2). 

The data were reduced by ESO personnel us- 
ing the FLAMES-UVES reduction pipehne (see 



measured RVs are presented in Table [T] with their errors; 
RVs for all observed stars are the weighted mean from the 
two measurements for each star, and these provide estimates 
of the errors by taken from f xcor. 



3 MEMBERSHIP AND CLUSTER 
KINEMATICS 

The derived radial velocities and met allicit ies, along with 
the available photometry of To2 st ars l|Kub iak ct al. I99J; 
IPhelps. Janes fc Montgomervl[l994l ). are used to revise the 
fundamental parameters of the cluster. Cluster membership 
for the observed stars was determined on the basis of their 
radial velocity (Table [l]). The heliocentric velocities for all 
the observed s tars (Vr) are report ed in Table [T] which pro- 
vides IDs from lPhelps et al.l l|l994 col. 1) and lKubiak et al.1 
(Too? ml 9^ stellar coordinates, magnitudes and colours 



http://www.eso.org/project s/dfs/dfs-shared/w eb/vlt/vlt-instrum Birt^]pfi^ft(gi^ti}V ^994) ^y^^^^ FigureE we can eas- 



for documentation on the FLAMES-UVES pipeline and 
software) which corrects the spectra for the detector bias 
and fiat-field. Then a wavelength calibration based on 
Th-Ar calibration-lamp spectra was applied. Finally, each 
spectrum was flux-calibrated by applying the response-curve 
of the instrument, and the echelle orders were combined to 
obtain a single mono-dimensional spectrum. The resulting 
UVES spectra have a dispersion of 0.1 A pixel" ^ and a 
typical S/N ~ 15-20. 

The GIRAFFE spectra were obtained simultaneously 
using the HR09B set-up {R = 21,000 in the 5138-5350 
A range) with 35 stars targeted. These data were re- 
duced similarly to the UVES spectra, are at a dispersion 
of O.25A pixer\ and have S/N ~ 60-70. 

The stars targeted with GIRAFFE and UVES were ob- 
served twice, so that two independent RV measurements 
were derived for each star from the collected spectra. Radial 
velocities were measured using the IRAF utility fxcor; this 
routine cross-correlates the observed spectrum with a tem- 
plate having known radial velocity. As a template, we used a 
synthetic spectrum calculated for a typical solar metallicity 
giant star [%s - 5000 K, log{g) = 3.0, vt = 1.3 km s"^]. 
Then each fxcor measured RV was converted to a helio- 
centric velocity using the IRAF routine rvcorr. All of the 



ily identify cluster stars from their RV. Stars with radial 
velocities of 121 ± 4 km s~^ are considered to be cluster 
members on the basis that: (1) examination of the RV distri- 
bution in Figure [T] clearly shows a peak associated with the 
cluster at 121 km s~^, (2) open cluster velocity dispersions 
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(2002). Thus, the ± 4 km s~^ acceptance window for mem- 
bers is more than generous for including suspected RV mem- 
bers, but still small enough to exclude most (but possibly 
not all) expected Galactic contamination. From the distri- 
bution of RVs shown in Figure 1, one might expect field 
star contamination across an 8 km s~^ window at the level 
of no more than a f ew stars. Moreover, from the Besancon 
model l|Robin|[2003l ) we expect the mean velocity of Milky 
Way stars along the To2 line of sight to be about 57 km s~^, 
with a dispersion of 38 km s~^; which suggests, given our 
sampling, a contamination rate of 0.02 stars (~ 2%) with 
RVs of 117-125 km s~^ within the observed field of view. 

Analyzing the RV membership sample (shown within 
the dotted line s in Figure [5} using techniques from 
iPrvor fc Mevlar] l|l993l ). we derive a systemic RV for To2 
of (Vr) = 121.04 ± 0.43 km s~^ and a total velocity disper- 
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Table 1. Radial velocities, coordinates, photometry, cluster membership, and atmospheric parameters for 
stars with sufficient S/N to derive abundances from the GIRAFFE and UVES stellar spectra. 
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sion of ai„t = 1.81 ± 0.30 km s'V Of course, this dispersion 
may be somewhat inflated by any field star contamination 
and/or motions of unresolved binary stars. Cluster member- 
ship and the To2 velocity dispersion is addressed further in 

m 

4 ABUNDANCE ANALYSIS 

4.1 Atomic parameters and equivalent widths 

The analysis of chemical abundances was carried out with 
the l atest version (2005) of the program MOQG tSnedeni 
1 19731 ) and using model atmospheres by iKurucj (| 19791 ). 
MOOG performs its analysis in local thermodynamic equi- 
librium (LTE). We derived equivalent widths of unblended 
spectral lines by using a semi-automatic SuperMongo pro- 



gram written by one of the authors. Repeated measure- 
ments show a typical error of about 5-lOmA for the weak- 
est lines because of the relatively low S/N of the spectra 
(UVES: ~ 20, GIRAFFE: 60-70 ) compared to those typi- 
cally used i n chemical abu n dance studies. The line list was 
taken from iGratton et all (|2003l ) for UVES spectra, and 
from the VALEQ (Vienna Atomic Line Database; Kupka et 
al. 1999) database for GIRAFFE spectra. The log(g/) pa- 
rameters of these lines were redetermined by a solar-inverse 
analysis measuring the equivalent widths from the NOAO 
solar spectrum and adopting the standard solar parameters 
[Teff = 5777 K, log(5) = 4.44, and Vt = 0.8 km s"^]. The so- 
lar abundances found by this analysis are reported in Table 
O For each element we give also the number of measured 

^ Available at |http://www.astro.uu.se/~vald/| 
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Figure 2. (a) The distribution of all new and previously measured radial velocities (with error bars) for stars in the field of Tombaugh 2 
as a function of their projected distance from the cluster centre. Squares denote stars observed in this study, with filled squares denoting 
those observed with GIRAFFE and open squ ares represen t ing th ose with UVES spectr a. Stars selected to be RV members are shown 
by black symbols. Triangles denote RVs from [Brown et ah! ll 19961 ) and iFriel et al.l ll2002h . with darker triangles denoting members from 
these studies by our RV definition, (b) Same as (a) but with metallicity information for member stars shown by colour: red points are 
"metal-rich" ([Fe/H] > —0.15) and blue points are "metal-poor" {[Fe/H] < —0.15) — see the discussion in 



Table 2. Adopted solar abundances for UVES and GIRAFFE 
spectra. 



Element UVES # lines GIRAFFE # lines 



Fe I 


7.48 


50 


7.51 


40 


Fe II 


7.51 


4 


7.52 


6 


Na I 


6.31(LTE) 


2 






Mg I 


7.53 


1 






Si I 


7.61 


1 






Ca I 


6.37 


2 


6.55 


2 


Ti I 


4.93 


9 


5.08 


7 


Ti II 


4.96 


11 


5.08 


3 


Cr I 


5.65 


5 


5.66 


7 


Cr II 


5.72 


2 






Ni I 


6.26 


11 


6.51 


4 


Ba II 


2.45 


2 







Table 3. Sensitivity of derived abundances to the atmospheric 
parameters. 



A(T) A(log(g)) Aivt) 

Element (per +100 K) (per +0.2) (per +0.2 km s^^) 



Fe I 


+0.07 


+0.02 


-0.07 


Fe II 


-0.07 


+0.10 


-0.06 


Na I 


+0.08 


-0.04 


+0.04 


Mg I 


-0.08 


-0.04 


+0.06 


Si I 


-0.02 


+0.02 


+0.05 


Ca I 


+0.09 


-0.04 


-0.09 


Ti I 


+0.15 


-0.02 


-0.11 


Ti II 


-0.01 


+0.10 


+0.12 


Cr I 


+0.13 


-0.05 


-0.10 


Cr II 


-0.04 


+0.09 


-0.10 


Ni I 


-0.08 


+0.00 


+0.11 


Ba II 


+0.03 


+0.05 


-0.15 



4.2 Atmospheric parameters 



spectral lines. In one case (Ni i) the solar abundances for 
UVES and GIRAFFE spectra diflter by more than 0.2 dex. 
The reason is that we decided to adjust the \og{gf) values 
only to remove the scatter affecting each line list and not to 
register the two line lists onto a common scale; as a result 
there remains a possible systematic offset for Ni i. 



Initial estimates of the atmospheric parameter T^a were ob- 
taine d from photometric observations using the relations 



Alonso et all (^1999^ . We first adopted E{B - V) = 
Brown et aLlll996l ) to correct colours for the interstel- 



from 
0.35 

lar extinction. We then adjusted the effective temperature 
to minimize the slope of the abundances obtained from Fe 
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I lines with respect to the excitation potential in the curve 
of growth analysis. Initial guesses for the gravity, log g were 
derived from the canonical formula: 



Table 4. Measured Abundances for UVES stars. 



log ^ 



log 



M 



+ 4 log 



rcff 



log 



L 



(1) 



In this equation the mass M/Mq was derived from 
the comparison between the position of the star in 
the colour-magnitud e diagram and the Padova isochrones 
l|Girardi et al.l 2000l ). The luminosity L/Lq was derived 
from the absolute magnitude My, assuming a distance mod- 
ulus of (m - M)v = 15.3 (Kubiak et al. 1992). The bolo- 
metric co rrection (BC) was d erived from the BC-Te^ rela- 
tion from lAlonso et al] (119991 ). The input logg values were 
then adjusted in order to satisfy the ionization equilib- 
rium of Fe I and Fe II during the abundance analysis. Fi- 
nally, the microturbulence velocity is given by the relation 
l|Houdasheh. BeU fc Sweigartll2000l ): 



vt = 2.22 - 0.322 log 5 



(2) 



We then adjusted the microturbulence velocity by minimiz- 
ing the slope of the abundances obtained from Fe I lines 
with respect to the equivalent width in the curve of growth 
analysis. From the derived spectroscopic temperatures, we 
are a ble to obtain the int rinsic {B — V) colours for our stars 
using IaIouso et ah! (|l999l ). We find that for these stars that 
E{B — V) ^ 0.30 ± 0.04, in agreement with the reddening 
estimate for To2 bv lBrown et al.1 (Il996l ). The final derived 
stellar parameters are listed in Table [1] The uncertainties 
in the spectroscopically determined stellar parameters are 
of the order of ±100 K in T^ff, ±0.2 in log{g), and ±0.2 
km s~^ in micr oturbolent velocity for red giant stars (see 
iFriel et al]|2005^ . Table [3] lists the impact of uncertainties in 
the atmospheric parameters on the derived abundances for 
the elements considered in our analysis. Variations in param- 
eters of the model atmospheres were obtained by changing 
each of the parameters one at a time for one of the analyzed 
stars (#179), assumed to be representative of all the stars 
considered in this paper. 



4.3 Stellar abundances 

The derived abundances from the UVES spectra are pre- 
sented in Table |4] In addition to the UVES spectra, we also 
derive some abundances from the bright GIRAFFE spec- 
tra of cluster members (membership is described above in 
§3), using the same techniques described above. The derived 
abundances from the GIRAFFE spectra are listed in Tables 
[S] and [51 together with their uncertainties. 

The Na abundance was obtained from the spectral lines 
at 5662-8 and 6154-60 A. These features are well known to 
be affected by NLTE effect s. For this reason we applied an 
NLTE correction from .Gratton et al.l (|l999l ) to the output 
LTE abundances (see Tabled}. 



5 CHEMICAL ABUNDANCES AND CLUSTER 
KINEMATICS REVISITED 

Our abundance analysis of To2 finds a surprising spread in 
[Fe/H] for stars selected as members based on their RVs. In 
Figure |3l we show the distribution of [Ti i/Fe] versus [Fe/H] 



ID 


El 


(X) 




EL0 


# 


[El/X] 


158 


Fe I 


(H) 


7.42±0.04 


7.48 


31 


-0.06±0.04 


158 


Fe II 


(H) 


7.45±0.14 


7.51 


3 


-0.06±0.14 


158 


Na 1(568) 


(Fc) 


6.51±0.07 


6.31 


2 


-f0.26±0.08^ 
-l-0.31±0.08'^ 


158 


Nai(616) 


(Fc) 


6.17 


6.31 


1 


-0.08-^ 
+0.07^ 


158 


Mg I 


(Fe) 


7.25 


7.53 


1 


-0.22 


158 


Si I 


(Fc) 


7.72±0.10 


7.61 


2 


+0.17±0.11 


158 


Ca II 


(Fe) 


6.34±0.07 


6.37 


5 


-f0.03±0.08 


158 


Til 


(Fe) 


5.07±0.07 


4.93 


7 


-f-0.20±0.08 


158 


Ti II 


(Fe) 


4.93±0.10 


4.96 


2 


+0.03±0.11 


158 


Cri 


(Fe) 


5.60±0.13 


5.65 


5 


-f0.01±0.14 


158 


Cm 


(Fe) 


5.64 


5.72 


1 


-0.02 


158 


Nil 


(Fe) 


6.21±0.08 


6.26 


7 


-l-0.01±0.09 


158 


Bai 


(Fe) 


1.98±0.08 


2.45 


2 


-0.41±0.09 


158 


a 


(Fe) 








-l-0.04±0.08 


162 


Fe I 


(H) 


7.44±0.04 


7.48 


35 


-0.04±0.04 


162 


Fe II 


(H) 


7.47±0.10 


7.51 


4 


-0.04±0.10 


162 


Nai(568) 


(Fc) 


6.87±0.09 


6.31 


2 


-l-0.60±0.10^ 
-f0.65±0.10^ 


162 


Nai(616) 


(Fc) 


6.35 


6.31 


1 


-fO.08-^ 
-1-0.23^ 


162 


Mgi 


(Fc) 


7.55 


7.53 


1 


-1-0.06 


162 


Si I 


(Fc) 


7.72±0.08 


7.61 


2 


-f-0.15±0.09 


162 


Ca II 


(Fc) 


6.26±0.09 


6.37 


5 


-0.07±0.10 


162 


Til 


(Fe) 


4.99±0.09 


4.93 


8 


-f0.10±0.10 


162 


Ti II 


(Fe) 


5.15±0.16 


4.96 


2 


-f0.23±0.16 


162 


Cri 


(Fe) 


5.80±0.13 


5.65 


3 


-f0.19±0.14 


162 


Cm 


(Fe) 


5.78±0.06 


5.72 


2 


-l-0.10±0.07 


162 


Nil 


(Fc) 


6.29±0.07 


6.26 


7 


-l-0.07±0.08 


162 


Bai 


(Fc) 


2.10±0.10 


2.45 


2 


-0.31±0.11 


162 


a 


(Fc) 








-l-0.09±0.05 


164 


Fe I 


(H) 


7.39±0.03 


7.48 


49 


-0.09±0.03 


164 


Fe II 


(H) 


7.42±0.10 


7.51 


3 


-0.09±0.10 


164 


Na 1(568) 


(Fc) 


6.36±0.07 


6.31 


2 


-f0.14±0.07^ 
-f0.19±0.08'^ 


164 


Nai(616) 


(Fc) 


6.23 


6.31 


1 


-1-0.01-^ 
-1-0.16^ 


164 


Mg I 


(Fc) 


7.30 


7.53 


1 


-0.14 


164 


Si I 


(Fc) 


7.62±0.05 


7.61 


2 


-f0.10±0.06 


164 


Ca II 


(Fc) 


6.30±0.05 


6.37 


9 


-f0.02±0.06 


164 


Til 


(Fc) 


5.14±0.04 


4.93 


11 


-l-0.30±0.05 


164 


Ti II 


(Fc) 


5.23±0.08 


4.96 


2 


-l-0.36±0.08 


164 


Cri 


(Fe) 


5.77±0.11 


5.65 


5 


-f0.21±0.11 


164 


Cm 


(Fe) 


5.89±0.07 


5.72 


2 


-f0.26±0.08 


164 


Nil 


(Fe) 


6.25±0.06 


6.26 


11 


-f0.08±0.07 


164 


Bai 


(Fe) 


2.52±0.01 


2.45 


2 


-l-0.16±0.02 


164 


a 


(Fc) 








-f0.13±0.10 


165 


Fe I 


(H) 


7.41±0.05 


7.48 


26 


-0.07±0.05 


165 


Fe II 


(H) 


7.44±0.09 


7.51 


4 


-0.07±0.09 


165 


Na 1(568) 


(Fc) 


6.36±0.01 


6.31 


2 


-f0.12±0.05^ 
-l-0.17±0.05'^ 


165 


Nai(616) 


(Fc) 


6.43 


6.31 


1 


-1-0.19-^ 
-1-0.34^ 


165 


Mgi 


(Fc) 


7.18 


7.53 


1 


-0.28 


165 


Si I 


(Fc) 


7.72±0.18 


7.61 


2 


-f0.18±0.19 


165 


Ca II 


(Fc) 


6.30±0.07 


6.37 


8 


-l-0.00±0.09 


165 


Til 


(Fc) 


5.18±0.11 


4.93 


7 


-l-0.32±0.12 


165 


Ti II 


(Fc) 


5.01±0.64 


4.96 


2 


-l-0.12±0.64 


165 


Cri 


(Fc) 


5.77±0.14 


5.65 


5 


-l-0.19±0.15 


165 


Cm 


(Fc) 




5.72 







165 


Nil 


(Fe) 


6.37±0.08 


6.26 


4 


-l-0.18±0.09 


165 


Bai 


(Fe) 


2.40±0.16 


2.45 


2 


-f0.02±0.17 


165 


a 


(Fe) 








-0.07±0.11 


^ NLTE solution 



L T.TP, snliitir.Ti 



Table 5. Measured Abundances for GIRAFFE stars. 



ID El (X) ELMea. EL0 # [El/X] 
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Table 6. Measured Abundances for GIRAFFE stars (Cent). 













4± 
W 




191 


Fe I 


(H) 


7.26±0.06 


7.51 


19 


-0.25±0.06 


191 


Fell 


(H) 


7.23±0.48 


7.52 


3 


-0.29±0.48 


191 


Til 


(Fe) 


5.21±0.18 


5.08 


4 


+0.38±0.19 


191 


Till 


(Fe) 


4.99±0.26 


5.08 


2 


+0.16±0.27 


191 


Cri 


(Fe) 


5.29±0.07 


5.66 


3 


-0.12±0.09 


191 


Nil 


(Fe) 


6.35±0.11 


6.51 


2 


+0.09±0.12 


196 


Fe I 


(H) 


7.23±0.06 


7.51 


23 


-0.28±0.06 


196 


Fe 11 


(H) 


7.26±0.08 


7.52 


3 


-0.26±0.08 


196 


Til 


(Fe) 


5.30±0.17 


5.08 


3 


+0.50±0.18 


196 


Ti II 


(Fe) 


4.93±0.16 


5.08 


3 


+0.13±0.17 


196 


Cri 


(Fe) 


5.03±0.25 


5.66 


4 


-0.35±0.26 


196 


Nil 


(Fe) 


6.31±0.36 


6.51 


4 


+0.08±0.36 


199 


Fe I 


(H) 


7.30±0.04 


7.51 


20 


-0.21±0.04 


199 


Fe II 


(H) 


7.31±0.42 


7.52 


3 


-0.21±0.42 


199 


Til 


(Fe) 


4.99±0.15 


5.08 


5 


+0.12±0.15 


199 


Ti II 


(Fe) 


5.23±0.06 


5.08 


2 


+0.36±0.07 


199 


Cri 


(Fe) 


5.33±0.15 


5.66 


4 


-0.12±0.15 


199 


Nil 


(Fe) 


6.50 


6.51 


1 


+0.20 


231 


Fe I 


(H) 


7.42±0.05 


7.51 


23 


-0.09±0.05 


231 


Fe II 


(H) 


7.43±0.11 


7.52 


2 


-0.09±0.11 


231 


Til 


(Fe) 


5.46±0.12 


5.08 


5 


+0.47±0.13 


231 


Till 


(Fe) 


5.11±0.10 


5.08 


2 


+0.12±0.11 


231 


Cri 


(Fe) 


5.32±0.11 


5.66 


5 


-0.25±0.12 


231 


Nil 


(Fe) 


6.35±0.10 


6.51 


3 


-0.07±0.11 




-0.4 -0.2 0.2 

[Fe/H] 

Figure 3. Plot of the derived [Fe/H] vs. [Ti l/Fe] for the "mem- 
bers" from our analysis. Red points denote "metal-rich" and blue 
points are "metal-poor", same as Fig lb.. The dotted line de- 
notes the "mean" [Fe/H] value for the MR sample (—0.06), the 
dashed line for the MP sample (—0.28), and the solid line for the 
combined sample (—0.14) for the cluster if all stars are members. 
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Figure 4. Spatial distribution of To2 stars from IPhelps et al.l 
1119941 ) photometry. Marked stars are "members" from our analy- 
sis, with red points denoting "metal-rich" and blue points "metal- 
poor." Open boxes denote the UVES stars. The gray circle de- 
notes a typical open clusters tidal radius of 10 pc llPiskunov et al.l 
l2008fl scaled to the distance of To2. However as To2 is in the outer 
disc and older, and therefore presumibly more massive than the 
typical open cluster, the shown tidal radius should be taken as a 
lower limit to the extent of the cluster. 
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Figure 5. Integrated counts of stars with RVs with 15 arcmin 
of the cluster centre. The dotted line represents all stars with 
RVs (Brown et al. 1996; Friel et al. 2000; ours), while the dashed 
line is all member stars. The black solid line represe nts all of the 
members found in our study (excluding the stars of lBrown et al.l 
1 19961 : iFriel et al.ll2002l ) . Finally, the red line show the distribution 
of our "metal-rich" members and the blue line our "metal-poor" 
members. 



abundance ratios, which shows a strong overall correlation 
but also an apparent divide at [Fe/H] = —0.15 dex. Even if 
the division into two populations is not real and the metal- 
licity variation is more continuous, dividing the sample into 
two parts still provides a handy way to explore properties 
as a function of metallicity. To investigate this distribution 
further, we select two sub-samples by using the apparent di- 
vide at [Fe/H] = —0.15 dex and denote these as the "metal- 
rich" group (MR; [Fe/H] > —0.15 dex; red points) and the 
"metal-poor" group (MP; [Fe/H] < —0.15 dex; blue points). 

The spatial distribution of the two sub-samples is shown 
in Figure U as well as in Figure [^Ja. We find that the MP 
group is more centrally concentrated, as shown by the in- 
tegrated counts in Figure [51 though both groups could rea- 
sonably belong to the cluster. Even if we account for the 
selection effect of how stars were targeted for spectroscopy 
(the dotted line in Fig. 5) the MP stars still show a sig- 
nificant central concentration. The summary of kinemat- 
ics and abundances for these sub-groups are listed in Ta- 
ble [H showing [Fe/H], [Ti i/Fe], [Cr/Fe], [Ni/Fe], as well as 
the mean velocity ((K-)) and velocity dispersion {(Jint). As 
may be seen from Table |8l the kinematical properties of the 
two samples are basically indistinguishable, so that kine- 
matics do little to help understand the origin of the metal- 
licity spread/differences among the RV-members. Figure [6] 
shows the distribution of MP and MR stars in the colour- 
mag nitude diagram (CM P) for To2, using the photometry 
from lPhelps et al] (|l994l ) and demonstrates that both pop- 



ulations also lie in areas of the CMD expected for member 
stars. 

5.1 Stellar abundance reliability 

Given that an abundance spread in an open cluster is an un- 
precedented find, one of the first points to explore is whether 
our derived metallicities, and the implied spread in chem- 
istry for stars with the RV of To2, could be artificially in- 
duced by errors in our analysis. Two possible sources of er- 
rors can mimic a spread in metallicity: 

(i) Differential reddening could alter colours an unknown 
amount, thereby affecting the derived temperature, produc- 
ing a spread in metallicity. However, for differential redden- 
ing to explain the observed abundance spread across both 
samples (0.4 dex) would require reddening differences of 
> 0.2 E{B — V), which is not possible given the observed 
CMD and our spectral analysis. Moreover, Figure 4 shows 
that the MP and MR stars are mixed in their spatial distri- 
bution. In any case for our study the temperatures derived 
from photometry were used only as first guess in an iterative 
procedure that allowed us to find the atmospheric parame- 
ters directly from spectral data, and this spectral analysis 
reveals a reddening spread of only 5E{B — V)^ 0.04 ( il4.2p . 
As a result, we can exclude differential reddening as a pos- 
sible source for the observed metallicity spread in To2. 

(ii) Random errors in temperature can introduce random 
errors in the derived metallicity. However, such errors are 
difficult to estimate for our To2 data both because of the 
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Table 7. Measured Abundance Ratios for UVES stars. 





tr from M4 


a from M4 


Element 


(GIRAFFE) 


(UVES) 


error (Fel) 


0.06 


0.05 


error (Fell) 


0.06 


0.05 


error (Til) 


0.06 


0.04 


error (Till) 


0.10 


0.06 


error (CrI) 


0.10 


0.05 


error (CrII) 




0.08 


error (Ni) 


0.07 


0.03 


error (Sil) 




0.05 


error (Cal) 




0.03 


error (Ball) 




0.03 



relatively small numbers of stars in our sample, and because 
of the suspected intrinsic metallicity spread of our stars. The 
determination of the random error from T^f / and that of the 
abundance determinations are coupled, so that if the metal- 
licity spread is true it would lead to an overestimation of the 
random error due to the T^f / uncertainty. To get around this 
problem, we therefore compare our derived T^f / uncertainty 
to another analyzed cluster without an abundance spread to 
evaluate the reliability of our abundance determinations. We 
estimated the whole error (the error due to the equivalent 
width measurement and given by MOOG -|- the error due 
to the Te// uncertainty) referring to a study on the globu- 
lar cluster M4 (Marino et al., in preparation). Some of the 
authors of the present study of To2 also derived chemical 
abundances for a large statistic sample of stars in M4 us- 
ing the same method adopted here on spectra having the 
same wavelength range and comparable S/N. Assuming a 
homogenous composition for the stars in M4 (an assumption 
usually true for most of the elements in a globular cluster 
with the exception of the light ones) a dispersion (a) was 
found in the iron content of 0.06 dex for GIRAFFE data 
and of 0.05 dex for UVES ones. These observed dispersions 
can be considered as a good estimate of the whole random 
error to be used in our study of To2. In Table [7] we report 
the dispersion for all the elements considered in this study 
as derived from the M4 data. Given that our measured un- 
certainties on the abundance determination are consistent 
between the M4 and To2 analysis, we find that errors in the 
determination of Te// cannot be the source of the observed 
abundance spread. 

We find that random observational errors cannot ex- 
plain the mean difference of 0.22 dex (see Table[S| in metal- 
licity observed in To2 between the metal rich and the metal 
poor group. Additionally, since systematic errors introduced 
by the method affect all the results in the same way, these 
also cannot explain the observed abundance spread. At this 
point we conclude that the observed spread in metallicity in 
our sample is real. 

5.2 Comparison to previous findings 

iBrown et al.l (|l996l ) conducted a high resolution (R ~ 
34,000) study of To2 with the CTIO 4-m telescope and found 
[Fe/H] = -0.4 ± 0.25 for E{B-V) = 0.4 or [Fe/H] = -0.5 ± 
0.25 for E{B-V) = 0.3 bas ed on three stars, wi th [Fe/H] = 
-0.2, -0.4, and -0.57. The lBrown et al.l l|l996l ) abundance 



analysis showed that To2 has a reddening of E{B — V) — 0.3- 
0.4, which is consistent with our spectroscopically measured 
reddeni ng. Unfortuna t ely w e do not have stars in common 
with the lBrown et al.| l| 19961 ) sample for a direct comparison 
of our metallicity derivations, but the metallicities of their 
three stars are within the range of our derived [Fe/H], ac- 
counting for the errors, for To2 members. However, the more 
metal-rich of the three Brown et al. stars is ruled out as a 
cluster member based on our radial velocity criterion. 

Another study of cluster metal licities including To2 
was conducted by iFriel et al.1 (|2002| ) using the CTIO 4- 
m/ARGUS, which yielded a much lo wer resolution (R ~ 
1300) than either our study or that of iBrown" et al.1 l| 19961) ■ 
and with metall icities determined from spectral indices. 
iFriel et al.l l|2002l ) found [Fe/H] = -0.44 ± 0.09 for To2 
from a sample of 12 member stars, and with t he individ- 
ual m easurements ranging from -0.28 to -0.65. iFriel" et al.l 
(|2002h found nearly identical meta llicities for the two stars in 
common w ith Brown et al.l (|l996l ). We observed three other 
iFriel et al.l (|2002l ) stars with the GIRAFFE spectrograph 
but find poor agreement with the previous low-resolution re- 
sults in 2 of the 3 st ars, as shown in Table[^ where I D is fr om 
iPhelps et all (|l994l ) and ID2 is from lKubiak et al.1 (Il992l ). It 
should be pointed out that sometimes high-resolution stud- 
ies to find slightly higher metallicities when compared with 
low-resolution studies based on spectral indices (e.g., in the 
high resolution study of Berkeley 17, Friel et al. 2005 find 
[Fe/H] = -0.10±0.09 as compared to [Fe/H] = -0.33±0.13 
from the low- resolution spectroscopy of Friel et al. 2002) . 



5.3 Revised reddening and distance 

From the stars that have been selected as members, we 
have determined the cluster mean metal abundance for the 
MP and MR sub-samples (see Table IHJ. We have generated 
isochrones for the exact metallicity of each cluster sample, 
transforming the mean [Fe/H] i nto Z, using Padova mod- 
els and following ICarraro et al. I 11999). The corresponding 
by-eye fit isochrone is then superposed on the CMD, with 
stars of the other abundance sample removed, and fitted 
to the data distribution (Figure [B]). In Figure [71 we match 
Padova isochrones to the {V, V—I) CMDs from lPhelps et al] 
(1994). We show these matches with stars from the inner 2 
arcmin of the clusters plus RV members selected from the 
MP (Figure [XK) and MR (Figure [7)d) samples. Estimates of 
the basic parameters (age, distance and reddening) are de- 
rived from comparison to the Padova isochrones. Our best 
fit for both samples is: age — 2.0 Gyr, E{B — V)— 0.3, and 
(m - M)o = 14.5 or d = 7.9 kpc and Roc = 14.2 kpc. Ad- 
ditionally, the spectroscopy yields a reddening in the range 
0.3-0.35 which is consistent with isochrone fitting to both 
samples; how ever, we find th e value of E{B — V) = 0.40 
found by Ku biak et al.l ([l99^) cannot provide a g ood fit to 
the (V, V-I) CMDs using the lPhelps et al.1 (|l994( ) photom- 
etry. Our spectroscopically determined red dening is also con- 
sisten t with a newer near-infrared study (|Kveong fc BvunI 
l2000h who found E{B ~ V) = 0.24 ± 0.12. 




Figure 6. (a) Colour-magnitude diagram for stars within a 2' radius of To2 from iPhelps et al. plus all stars with current and 

previous RV determinations. Squares d enote observations with VLT/FLAM ES+GIRAFFE, open squares with VLT/FLAMES+UVES, 
and triangles are previous observations llBrown et al ] |l996l : iFriel et al.|[20o3) . Red and blue symbols arc members based on their RV. All 
symbols same as Figure 1. (b) RGB and red clump from (a) with a Padova isochrone of Z = 0.0164, age = 2.0 Gyr, d = 8.4 kpc, and 
E{B -V) = 0.25 overplotted. 



Table 8. GIRAFFE Spectra Radial Velocity and Abundance "Sub-Populations" 



Sample 


# Stars 


[Fe/H] 


[Ti l/Fe] 


[Cr/Fe] 


[Ni/Fe] 


(Vr) 
















(km s ^) 


(km s"-*-) 


Metal-rich (UVES only) 


4 


-0.07±0.01 


-l-0.23±0.06 


-l-0.15±0.05 


-l-0.08±0.04 


119.2±1.0 


2.1±0.7 


Metal-rich (GIRAFFE only) 


7 


-0.05±0.02 


-l-0.14±0.05 


-0.48±0.12 


-0.20±0.09 


121.7±0.5 


1.2±0.3 


Metal-rich (UVES-I-GIRAFFE) 


11 


--0.06±0.01 


-f0.23±0.03 


-0.25±0.07 


-0.05±0.05 


120.8±0.6 


2.0±0.4 


Metal-poor (GIRAFFE) 


7 


-0.28±0.03 


-l-0.32±0.06 


-0.26±0.06 


-0.01±0.07 


121.5±0.5 


1.4±0.4 


All RV Members 


18 


-0.14±0.03 


+0.24±0.03 






121.0±0.4 


1.8±0.3 



6 DISCUSSION 

While this study was intended to verify or refute the connec- 
tion of To2 with the Monoceros stream/Galactic anticenter 
stellar structure (GASS), our work has produced more ques- 
tions than answers. The observed spread in To2 metallicity 
and Q-abundances was not expected and cannot be simply 
explained. To date, no other open cluster is known to have 
a metallicity spread, and it is difficult to understand how 
a low-mass open cluster might have been able to retain gas 
after an initial starburst to self-enrich. Next, we review our 
findings and possible scenarios that could explain the data 



6.1 Are both populations part of Tombaugh 2? 

First, we ask the question of whether both the MP and MR 
groups among the To2 RV members are likely to be truly 
part of the cluster, and then investigate the likelihood that 



either population is simply contamination of the RV member 
sample. 



6.1.1 Is the metal poor population part of Tombaugh 2? 

To2 has been prop osed to be part of the putative 
GASS cluster system l|Frinchabov et al.|[20oi : iMartin et all 
l2004al '). and the MP population is similar, though slightly 
more metal rich than, other studied Monoce ros candi- 
date clusters such as Saurer 1 (iFrinchabov fc Phe lps 200^; 



ICarraro fc BaW '2003': [C arraro et al.l 



Carraro et all .2004: ,Yong et all |2005|) 
Rosenberg et al.l 19981 : ISaraiedini et al 



2004D . Berkeley 29 
and Palomar 1 
|2003). These clus- 
ters all have ages of approximately 4-6 Gyr, metallicities 
between —0.7 and —0.4 dex, and are all Q-enhanced. Our 
To2 MP group is younger though similarly moderately metal 
poor ([Fe/H] ~ —0.3) and also a-enhanced. This is at least 
circumstantial evidence that the MP group might be the 
more likely member population. This subgroup also has a 
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Table 9. Comparison of lFriel et all l|2002l ) and GIRAFFE Data 



ID 


ID2 


yr,GIRAFFE 




A Vr 


[Fe/H]oiB.AFFE 


[Fc/R]Fr^el 


A[Fe/H] 






(km S~^) 


(km s"-*-) 


(km s"-*-) 








63 


20 


120.96 ± 0.16 


123 ± 10 


2 


-0.19±0.05 


-0.40±0.10 


0.21 


98 


28 


121.71 ± 0.18 


132 ± 10 


10 


-0.06±0.05 


-0.48±0.22 


0.42 


177 


61 


118.74 ± 0.16 


139 ± 10 


20 


-0.43±0.05 


-0.39±0.31 


0.04 




Figure 7. CMDs for To2 showing best fit Padova isochrones. (a) Padova Z = 0.00955 ([Fe/H] = -0.28) isochrone fit is shown with MP 
members and other RV members without metaUicity determinations, (b) Padova Z = 0.01665 ([Fe/H] = —0.06) isochrone fit is shown 
with MR, members and other RV members without metaUicity determinations. 



metaUicity roughly cons istent with p r eviou s met aUicity de- 
termin ations for To2 bv lBrown et al.l (Il996l ) and lFriel et al.l 
— though it is not impossible that these other re- 
searchers could have also been hit with a spate of bad luck 
in terms of misleading target selection. 

On the other hand, it must also be admitted that the 
MP group both exhibits more scatter in terms of metaUic- 
ity properties than the MR population and has a smaller 
number of members — though this number increases when 
we account for the addi tional "MP " mem bers found in the 
iBrown et all l| 19961 ) and lFriel et~aLl (|2002l) samples that are 
not part of our sample. Perhaps the strongest support for the 
supposition that the MP group is a "real" part of the cluster 
comes from the fact that the MP group is strongly centrally 
concentrated (Figs. 4 and 5) and that it has a marginally 
smaller velocity dispersion than the MR group, as shown in 
Table H 

6.1.2 Is the metal rich population part of Tombaugh 2? 

The MR population is indeed less centrally concentrated 
that the MP population, and in fact has a distribution 
function not unlike our spectroscopic sampling function, as 
shown in Figure |4] and [5] Given the small number statistics 



overall, it may not be impossible for us to have been success- 
ful at achieving a uniform "success" fraction (i.e. in iden- 
tifying RV members) over the radial range sampled, even 
by chance. Indeed, the overall distribution of RV members 
very closely tracks the sampling function. Though Figure |4] 
clearly shows a tighter MP distribution than MR popula- 
tion, it is not inconceivable for there to exist a metaUicity 
gradient in the system to yield a differential radial distribu- 
tion — on the other hand, the one needed (i.e. more metal 
poor towards the centre) is quite unusual and not seen else- 
where. 

On the other hand, the MR group has more members 
and is more tightly clumped in [Fe/H]. Given the tight 
RV distribution, chemical properties, and sheer numbers it 
seems difficult to believe that the MR population is not some 
coherent Simple Stellar Population (SSP), whether part of 
the cluster or a "contaminant." We show below that this 
population is quite unlike that of the disc field star popula- 
tion (which is not an SSP). 

If the true To2 stars are those in the MR group (and 
excluding the MP group) then To2 is less consistent with 
being part of the putative GASS cluster system, as it would 
be inconsistent with the observed metaUicity distribution of 
Monoceros (Crane et al 2003; Sbordone et al 2005; de Jong 
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et al. 2007). In this case, To2 would be the most distant 
metal-rich Milky Way disk open cluster found to date. It 
seems unlikely that a 2 Gyr old cluster with near solar [Fe/H] 
could have formed "normally" in the outer disc, given the 
observed, more depleted chemistry of the surro unding disc 
stars fe.g.. lCarnev et aLllioOsI : lYong et al.ll2006l ): this would 
suggest that the cluster must have formed elsewhere (e.g., in 
a dwarf galaxy) and been deposited in the outer disc. This is 
not out of line with current CDM models for the formation of 
Milky Way-like galaxies, which suggest that discs form from 
the outside via continued merging (e.g., Abadi et al 2003). 
Stars as metal rich as our MR group are indeed found in 
paradigm accretion events, like the Sagittarius dwarf (e.g., 
Smecker-Hane & McWilliam 2002, Monaco et al. 2005, Chou 
et al. 2007). 



6.1.3 Can either population simply he contamination? 

If we believe that only one of the two metallicity populations 
is truly part of To2, then the other population must be a 
contaminant of some kind. Here we address how likely it is 
that either population may be part of the Galactic disc. 

We have already seen (JJ)) that the predicted contam- 
ination rate of normal disc stars into our RV membership 
range based on the Besancon model (see Q is already 
smaller (by factor of 10^) than either our MP or MR sam- 
ples — thus, it seems very unlikely that either population 
can be a statistical, "fluke" superposition of a large number 
of disc stars with just the right RV. From Poisson statistics, 
the chances that we would get 7 (e.g., MP) or 11 (e.g., MR) 
stars as contaminants when the models predict 0.02 stars 
is 1.1 X 10~^ % and 2.3 x 10"" %, respectively. And this 
ignores the fact that other studies have identified additional 
stars in the To2 that we would categorize as "MP." 

Based on the radial distribution of members, the MP 
population seems more likely to be truly part of the clus- 
ter. This could then imply that the MR group is a con- 
taminant population. However, stars this metal rich (nearly 
solar metallicity) at Rgc > 15 kpc are not consistent with 
the measured metallicities of red giants in the outer disc 
Carney et al] l2005h as well as outer disc Cepheid stars 



Yong et al.ll20od ). which both suggest that the median disc 



metallicity at this Galactocentric distance should be [Fe/H] 
~ -0.4. These studies of outer disc tracers find no stars as 
metal rich as our MR group, even for the younger Cepheid 
populations. 

This then might suggest that the MP group could be a 
contaminating field population, since the chemistry of this 
gro up is more in line with that of the ICarnev et al.l (|2005l ) 
and lYong et al.l (|2006l ) stars in the outer disc. In turn, this 
implies that the tight radial concentration of the MP stars 
is a statistical anomaly. However Monte Carlo simulations 
show that the probability for a random disc population to 
achieve the observed level of concentration, taking into ac- 
count our sampling function, is 0.1%. It also implies that 
previous spectroscopic studies had similarly "bad luck" in 
picking targets to represent To2. 



6.2 Does Tombaugh 2 really have a metallicity 
spread? 

The cumulative evidence for To2 — the nearly identical ve- 
locities and velocity dispersions of the stars when divided 
into MP and MR groups, the positions of the stars in the 
CMD, the unlikeliness that these stars are field contamina- 
tion, etc. — suggests that all of the RV "members" in our 
To2 spectroscopic sample could truly be members of the sys- 
tem. If all RV "members" are indeed members of To2 then 
we are finding evidence for a significant metallicity spread 
within To2. That we also seem to find a similar bifurcation 
of [Ti i/Fe] with [Fe/H], rather than a scatter-diagram in 
Figure 3 is further evidence of a possible coherent chemical 
enrichment history. It is also possible, based on an apparent 
"break" in both [Fe/H] and [Ti i/Fe], that the To2 chemical 
and age distribution suggests the possibility that there may 
be two populations rather than a "trend." 

If these populations of varying chemistry do belong to 
To2, what formation scenarios could create such a unique 
open cluster? 



6.2.1 Is To2 ac 



two "very close" clusters? 



What is the possibility that two clusters with the same RV 
lie at nearly the same distance and have nearly the same 
metallicity and could also be aligned along the same line 
of sight? Though such a configuration could produce the 
observed results in the To2 field, and even account for the 
apparent "bimodal" behavior of the chemical distribution, it 
seems extremely unlikely given the small number of clusters 
(less than 20) known this far from the Galactic centre. It 
is the case that a superposition of clusters would increase 
the "discoverability" of this distant "double cluster", but 
the individual clusters would have to be very close together 
along the line of sight given that a significant broadening of 
the main sequence is not seen in the available photometry. 
However new higher quality photometry would be needed 
to fully test this case (as was don e with the Hubble Space 
Telescope in the case of NGC 2808: |Piotto et al.ll2007h . 

There are a few cases of overlapping clusters, with the 
most-well studied overlapping pair being NGC 1750 and 
NGC 1 758, with ages of 200 ± 50 and 400 ± 100 Myr respec- 
tively (iGaladf-Enn'quez et al.l Il998bl ). For NGC 1750 and 
NGC 1758 proper motions were required to be able to sep- 
arate the two clusters separated by ~ 130 pc along the line 
of sight (|Galadi-Enrfauez et al.iri998lj ). but with little no- 
ticeable differences seen in the overlapping CMD sequences 
(|Galadf-Enrfauez et al.|[l998al ). These clusters are the old- 
est pair of confirmed overla pping clusters, and, given the 
age and distance differences, iGaladf-Enriguez et al.l l|l998bl ) 
conclude that the clusters do not constitute a binary sys- 
tem. However for To2, the fact that mean velocity derived 
from high-resolution RVs is the same, and different than the 
Galactic disc trend, makes the possibility of it being part of 
a "double cluster" remote. 



6.2.2 Does To2 represent the merging of two clusters? 

Another explanation for our results in the To2 field might 
be that two clusters with nearly the same age and different 
metallicities could have have collided and merged. Given 
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that open clusters rotate with the Galactic disc, the small 
relative velocity difference between two clusters with similar 
orbits makes prolonged near encounters between close clus- 
ters possible. Additionally, the relatively large velocity dis- 
persion of To2, aint = 1-8 km s~^ (e.g., M67 has aint ~ 0.8- 
0.96 km s~^; Girard et al. 1989; Frinchaboy & Majewski 
2007) might suggest that the cluster could have been dynam- 
ically heated at some point, possibly by a merger of clusters. 
A simple calcula tion based on merging of two equal-mass 
eliptical galaxies (iBinnev fc Tremainelll987n . shows that for 
a two-body encounter with two equally massed clusters (as- 
suming each is 10* ^q) finds that mergers are possible if 
two clusters pass within less than 2 pc given a differential 
velocity of 5 km s~^. However as pointed out in the argu- 
ment above, given the small number of distant clusters the 
chances of such an occurrence are very small; nevertheless, 
the merging of two clusters cannot be ruled out as a possible 
explanation for the observed chemical peculiarities. 

6.2.3 Has To2 had multiple episodes of star formation? 

The cumulative evidence for To2 — the nearly identical ve- 
locities and velocity dispersions of the stars when divided 
into MP and MR groups, the positions of the stars in the 
CMD, the unlikeliness that these stars are field contamina- 
tion, etc. — suggests that all of the RV "members" in our 
To2 spectroscopic sample could truly be members of the sys- 
tem. If all RV "members" are indeed members of To2 then 
we are finding evidence for a significant metallicity spread 
within To2. That we also seem to find a similar bifurcation 
of [Ti i/Fe] with [Fe/H], rather than a scatter-diagram in 
Figure 3 is further evidence of a possible coherent chemical 
enrichment history. It is also possible, based on an apparent 
"break" in both [Fe/H] and [Ti i/Fe], that the To2 chemical 
and age distribution suggests the possibility that there may 
be two populations rather than a "trend." 

No other open cluster has been found to have an inter- 
nal metallicity spread. On the other hand, a small number 
of globular clusters have been found with metallicity spreads 
(e.g., uj Gen; Norris, Freeman, & Mighell 1996; Majewski et 
al. 2000; Garraro & Lia 2000; Frinchaboy et al. 2002; Vil- 
lanovaet al. 2007), multiple intern al populations (NG G 1851 
and NGG 2808; Milone et al. 2007. |Piotto et al.ll2007l '). or to 
represent one of a series of populations in a more complex 
structure (M54 in the Sagittarius dwarf spheroidal galaxy; 
Sgr dSph; Siegel et al. 2007). Interestingly, such unusual 
clusters typically have been associated with the possibility 
of tidal accretion events. 

Gould To2 be the remains of something that was once 
larger — a globular cluster or dwarf galaxy? This seems 
like the only way that To2 could have had two episodes of 
star formation, since the mass of a typical open cluster is 
generally much too small to retain and self-enrich gas after 
an initial star burst. The mass of a long-lived open cluster 
is typically a fe w lO'^-lO'^ MfT^ (e.g. , M67 has a current mass 
- 2 X 10^ M0; iHurlev et atlbOOSh . A mass of at least 10^ 
Mq would be needed for To2 to have retained enough gas to 
have a second epoch of star formation (e.g., the current mass 
of a; Gen is 2.8 x lO''; D'Souza & Rix 2005). Thus, this would 
require To 2 either to have had a significant dark matter 
content, which is ruled out by the velocity dispersion (and in 
any case would have made To 2 unique among star clusters. 



open or globular), or to have been initially much larger and 
then stripped down to its current size. Fortunately for this 
scenario, the disc is a dynamically brutal environment for 
clusters and dwarf galaxies, which are expected to undergo 
vigorous stripping. Indeed, this is the prevailing scenario to 
explain the presently visible multiple populations in the uj 
Gen "globular cluster" system. 

The two apparent sub-populations in To2 do seem to 
show a chemical trend evoking self-enrichment, with the 
metal-poor, a-enhanced population formed first, followed by 
formation of a more metal-rich, non-a-enhanced population. 
The above "evolution" is what is expected from closed-box 
chemical evolution models. Unfortunately, two other obser- 
vational facts contradict the self enrichment scenario: (1) 
The derived isochrone ages for the two populations are the 
same, within the errors of the fitting, which are ± 0.5 Gyr. 
(2) It is typical in self- enrichment scenarios for the more 
metal rich population, formed in the bottom of the poten- 
tial well more recently, to be more centrally concentrated. 
Thus, if To2 was the tidal remnant of a larger stellar system 
one would expect that the "younger" metal-rich population 
would be more centrally concentrated, as is found in the 
case of u) Gen (Pancino et al. 2000) and dSph galaxies (e.g.. 
Sculptor; Tolstoy et al 2004). However given the available 
data, we find that the metal-poor stars are more centrally 
concentrated, in marked contrast to large systems undergo- 
ing tidal stripping (Tolstoy et al 2004). 

6.2.i A "hybrid" solution? 

We have seen how a self-enrichment scenario seems to be 
contradicted by the relative radial and age distributions of 
the MP and MR populations. And we have seen that it is un- 
likely that To2 represents either dynamically or apparently 
merged/overlapped systems. Yet the strongly matched ve- 
locities and velocity dispersions of the two systems suggests 
a dynamical connection of the two populations, and there 
is at least some metallicity and velocity mismatch of these 
populations with the surrounding disc stars. There remains 
one possible scenario that could account for most of the ob- 
served properties in the To2 field, and which has a known 
prototype: The Sgr dSph stream and its cluster system. 

Any small spectroscopic survey of a Sgr star cluster, 
particularly one of those Sgr clusters near the core of the 
dSph (e.g., Terzan 7, Terzan 8 or Arp 2), would produce 
a number of the same phenomena we observe in our small 
spectroscopic survey of the To2 field: A centrally concen- 
trated, older core of more metal-poor stars that are part of 
the star cluster, seen against a more broadly distributed, 
typically younger and more metal rich population of stars, 
generally unbound, tidally stripped from the parent dSph 
galaxy, and all with the same radial velocity. This could be 
an analogous situation to what we observe in the field of To2, 
with the primary difference being that we are explaining the 
properties of an open, not globular, cluster. 

Thus, we hypothesize this as the leading explanation for 
what we are finding in the To2 field: We propose that the 
MP population represents the true cluster, as most strongly 
supported by the radial concentration of these stars, and 
that the MR population represents stars that are part of a 
distinct parent system, probably a tidally disrupted dwarf 
galaxy, that are spread more uniformly across our survey 
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field (and beyond). Such a scenario is consistent witii pre- 
vailing models of disc formation, which indeed suggest that 
discs continually grow by accreting "subhalos" around the 
edge (e.g., Abadi et al 2003, Brook et al. 2004). It also ac- 
counts for how we can have two populations with identical 
RVs, but differing metallicities, and with the more metal- 
rich population being more extended. That To2 has been 
already hypoth esized to be a part of the GASS/Monoceros 
"tidal stream" (|Frinchabov et al.ll2003 : lMartin et al.ll2004al ') 
provides additional support for this proposed scenario, al- 
though the velocity dispersion of GASS/Mon is much larger 
than our MR population (Crane et al. 2003, Martin et al 
2006). Nevertheless, it is possible that other tidal streams 
could exist and have contributed star clusters to the outer 
disc. 



6.3 Final remarks 

While we have discussed a number of possible explanations 
for our unusual findings in our survey of the To2 field, the 
only hypothesis that seems to explain most of the results in 
a natural way is that To2 is represented by our MP popu- 
lation and that the MR population represents stars from a 
parent dwarf galaxy that is contributing To2 to our Galac- 
tic disc. All other hypotheses — self-enrichment, a merger 
of two clusters, a superposition of two clusters — have seri- 
ous problems or are very improbable. It is also yet unclear 
whether To2 may be part of the proposed GASS/Mon sys- 
tem. Surely it is worth mentioning that if To2 actually has 
two different populations, then it is truly an unusual object. 
The one clear finding from this work is that further work is 
needed on the outer Galactic cluster system, in particular 
To2. A much larger spectroscopic sample that includes both 
high precision velocities and detailed abundances would be 
especially useful. Moreover, the properties of outer disc stars 
are needed to better constrain the existence, characteristics 
and origin of the GASS/Monoceros system and other tidally 
accreted systems in the outer Galactic disc. While RVs and 
abundances may not provide an explicit separation of an ac- 
creted satellite from the outer disc, it is nonetheless essen- 
tial that the kinematical and more importantly abundance 
trends in the outer disc are known. Without clear knowl- 
edge of the abundance patterns with radius, especially at 
large distances and in the second and third Galactic quad- 
rants, one cannot explicitly use abundances to distinguish 
whether the GASS/Monoceros stream is distinct from the 
outer Galactic disc. The present paper is one among many 
efforts that aim to correct the deficiency of abundance in- 
formation known for the outer Milky Way disc. 
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